Fission Cycling in Supernova Nucleosynthesis: Active-Sterile Neutrino Oscillations 
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We investigate nucleosynthesis in the supernovae post-core bounce neutrino-driven wind envi- 
ronment in the presence of active-sterile neutrino transformation. We consider active-sterile neu- 
trino oscillations for a range of mixing parameters: vacuum mass-squared differences of 0.1 eV^ < 
Srn^ < 100 eV'^, and vacuum mixing angles of sin'^ 2S„ > 10"**. We find a consistent r-process 
pattern for a large range of mixing parameters that is in rough agreement with the halo star CS 
22892-052 abundances and the pattern shape is determined by fission cycling. We find that the 
allowed region for the formation of the r-process peaks overlaps the LSND and NSBL (3-1-1) allowed 
region. 

PACS numbers: 14.60Pg, 14.60St, 26.30-|-k 



I. INTRODUCTION 

Neutrinos dominate both the dynamics and the nu- 
cleosynthesis in the region near the protoneutron star 
during a core collapse supernova. The prodigious neu- 
trino flux generated here is thought to deposit enough 
energy to liberate material from the gravitational well 
of the protoneutron star. Neutrino capture produces a 
neutron-rich environment outside the protoneutron star 
at late times post-core bounce. Understanding the role 
of neutrinos in this environment directly leads to our un- 
derstanding of heavy element nucleosynthesis in the core 
collapse supernova environment. 

One of the synthesis processes that may take place in 
this environment is the rprocess, or rapid neutron cap- 
ture process, which generates about half of the elements 
with atomic weight A > 100 P]. At the heart of the 
rprocess is the rapid capture of neutrons, which occurs 
much faster than the competing beta decay, forming very 
neutron rich nuclides. As the supply of free neutrons 
is exhausted, these nuclides decay back to beta stabil- 
ity, forming a characteristic r-process abundance pattern 
that can be observed today. 

While this basic r-process mechanism is understood, 
finding the astrophysical location of the rprocess proves 
more elusive. Several observational factors point towards 
the neutron-rich material produced near the protoneu- 
tron star of a core-collapse supernovae as a likely can- 
didate, e.g 0, Q • Sneden and Cowan 2003 y\ con- 
cluded the r-process abundance patterns in extremely 
metal poor giant stars match the second and third peaks 
of the solar system r-process pattern, indicating these el- 
ements were formed early in the evolution of the universe 
given these considerations Q . The neutrino-driven wind 
of the core collaspe supernova is a promising candidate 
for the rprocess, for a review see Wanajo et al. 2001 
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The evolution for a mass element which will eventu- 
ally undergo an rprocess in the neutrino-driven wind of 
theprotoneutron star proceeds through several stages 
la, H 113 • Material first emerges from the surface of 
the protoneutron star as free nucleons, carried off by the 
neutrino-driven wind. The strong influence of neutrinos 
in this regime produces neutron rich material resulting 
from the equilibrium effects of the i>e and Ve fluxes, as 
the Fg's are thought to have a longer mean free path at 
the surface of the protoneutron star and hence are more 
energetic. As the mass element moves farther from the 
star, it reaches lower temperatures (T < 750 keV), the 
nucleons coalesce into a particles and a large abundance 
of free neutrons. The a particles then combine into seed 
nuclei for the rprocess with 50 < A < 100. Neutron cap- 
ture begins at even lower temperatures (T < 300 keV) 
allowing the formation of r-process elements. 

It has long been suggested that fission could influence 
the rprocess primarily for cosmochronometers e.g 
|l2lll3l ll4| In sufficiently neutron-rich conditions, heavy 
elements which are unstable to fission can be produced, 
terminating the r-process path. Additionly, fission in- 
fluences the rprocess through the subsequent cycling of 
these fission unstable heavy elements by neutron capture 
on the fission products. Beta-delayed, neutron induced, 
and spontaneous fission are all thought to play important 
roles in determining the outcome of the final r-process 
abundances [Tsj l. 

Early r-process models in the neutrino-driven wind en- 
vironment, such as Woosley, et al, 1994 0, accurately 
reproduced the r-process abundances without significant 
fission cycling; However, later models produced lower en- 
tropy per baryon, producing nearly but not quite all the 
requisite conditions. Further, proper consideration of the 
near inertness of alpha particles to neutrino interactions 
resulted in decreased neutronization, presenting a signifi- 
cant impediment to r-process production in this environ- 
ment 11 El m El Ea 113. 

There are three possible solutions for circumventing 
these problems. One solution is modifying the hydrody- 
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namics by using a very fast dynamical timescale [l^ I2H 
even in a proton-rich environment , or increasing tem- 
perature with high entropy 0; however, it is not known 
how neutrino heating could generate such conditions. 
Another possibility is choosing a different location for the 
rprocess to occur, such as a neutron-star merger |2^l2^. 
However this solution is not currently favored by analysis 
of the observational data Q. Other astrophysical sites 
such as gamma ray burst accretion disks are currently 
being considered [25ll26| . 

The solution examined here is the introduction of 
active-sterile neutrino oscillations through the I'e ^ Vg 
and Ve ^ Vg channels. Two different types of terrestrial 
experiments place constraints on v,, ^ and Ve ^ Vg 
transitions, reactor experiments and short baseline accel- 
erator experiments. Relevant data of the former comes 
from CHOOZ [23| with it's null result on Ve disappear- 
ance, and the latter KARMEN, [2i|, with a nuU result 
on V ^ Ve appearance, leading to no indication of oscil- 
lations. LSND 29] is the only short baseline experiment 
which reported a positive signal for Uf, — > and this 
data is used as motivation for a sterile neutrino. 

Big Bang Nucleosynthesis(BBN) and the cosmic mi- 
crowave background (CMB) can provide constraints on 
the existence of sterile neutrinos. The number of new 
neutrino flavors, which would include a sterile neutrino, 
is restricted by the actions of the relativistic energy den- 
sity limits during the BBN epoch while simultaneously 
accounting for the WMAP data, e.g. '^HH. This re- 
striction may be lifted through flexibility in the electron 
neutrino asymmetry, which allows room for the addi- 
tion of neutrinos beyond the standard model [s^- CMB 
also places bounds on the neutrino mass; however, these 
bounds are dependent on the thermal spectrum of the 
sterile neutrino [3^ . At present, while they provide re- 
strictions, the BBN and CMB analysis do not discount 
light sterile neutrinos. 

We organize the paper as follows: In Section II, we 
describe the neutrino driven wind model and general nu- 
cleosynthesis in the neutrino-driven wind. Section III 
describes matter-enhanced active-sterile neutrino oscilla- 
tions and their effects on nucleosynthesis. Section IV de- 
scribes the network calculation of element abundances. 
In Section V we discuss the influence of fission in the 
rprocess and how fission cycling links the second and 
third peaks of the rprocess. Our results appear in Sec- 
tion VI showing our calculation of nucleosynthesis with 
active-sterile neutrino oscillations and we compare with 
both the solar system and halo star abundances. We 
conclude in section VII. 



II. NUCLEOSYNTHESIS IN THE NEUTRINO 
DRIVEN WIND 

The neutrino-driven wind epoch begins several seconds 
after post-core bounce in the type II supernovae environ- 
ment. The newborn protoneutron star undergoes Kelvin- 



Helmholtz cooling, radiating neutrinos that can lift mate- 
rial off the star's surface. This material is likely ejected 
by neutrino heating through charged-current reactions. 
It takes about ten neutrino interactions per nucleon to 
eject material from the surface of the star, and the neu- 
trinos set the electron fraction, Yg [loj |. 

In our calculations we use a one-dimensional model 
to follow a mass element's peregrination from a nucleon 
at the surface of the protoneutron star to an ejected 
r-process element. Although multidimensional phenom- 
ena such as convection may operate in this environment, 
this simple model elucidates the effect of neutrino flavor 
transformation on nucleosynthesis. 

A parameterization of the neutrino-driven wind in the 
type II supernovae environment following from ITsl 
l34 35, j56, J2j 2^ results in a constant mass outflow rate 
{M = Airr^pv) with a homologous outflow with a radial 
velocity v of 

war, (1) 

where r is the radial distance from the center of the 
protoneutron star. From this we find 

r = roexp(t/r), (2) 

where r,, is the initial radius and the expansion time 
scale, T is defined as 

T = r/v. (3) 

Including oc 1/r scaling for the enthalpy per baryon, 
the density p scales as 

p ex poexp{-3t/T). (4) 

In our calculations wc primarily use r — 0.3 s and 
S/k — 100 unless otherwise stated. Our mass element 
leaves the surface of the protoneutron star in Nuclear 
Statistical Equilibrium (NSE) due to the very high tem- 
perature. At entropies of S/k = 100 baryonic content 
exists only as protons and neutrons (35|. Electron neu- 
trino and anti-neutrino capture on free nucleons set the 
electron fraction and at the surface produce neutron-rich 
conditions. As material flows away from the protoneu- 
tron star it cools and once material reaches temperatures 
below T < 750 keV the formation of a particles begins. 

During this epoch since neutrino interactions have cre- 
ated a neutron- rich environment, all the protons and 
many neutrons bind into a particles. Electron neutrinos 
capture on the remaining neutrons producing protons, 
which in turn consume more neutrons to form additional 
a particles. This drives the electron fraction towards 
Ye — 0.5 decreasing the number of free neutrons and cre- 
ating an environment unfavorable to the rprocess, known 
as the a-effect [Hi^l. 
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FIG. 1: In the left panel (Fig. la) we show the contours of electron fraction, Ye, at the onset of rapid neutron capture, 
Tkev ~ 200. For the astrophysical conditions below, one expects an rprocess for Ye < 0.35 and an rprocess which matches the 
observed ratio of peaks heights for Ye < 0.3. For comparison with previous work, in the right panel (Fig. lb), we show the Ye, 
earlier, at the onset of heavy element formation, Tkev ~ 600. Neutrino captures on free nucleons increase the electron fraction 
slightly, between T^ev ~ 600 and Tkev ~ 200. The neutrino-driven wind parameters are entropy per baryon of S/k = 100, 
expansion time scale of r = 0.3 s, neutrino luminosities of — 10^^ ergs and La = 1.3 x 10^^ ergs, and neutrino temperatures 
T„ = 3.5 MeV and Ti, = 4.5 MeV. 



As the mass element continues moving away from the 
protoneutron star, the lower energies allow a particles 
and neutrons to combine into seed nuclei for the rprocess. 
Farther away from the protoneutron star, as the material 
falls out of NSE, it undergoes a series of quasi-equilibrium 
phases [13, |33, 113 ■ The large Coulomb barriers eventu- 
ally become insurmountable and charged current reac- 
tions drop out of equilibrium, while a large abundance 
of free neutrons remains. The mass element continues to 
move outward leaving neutron capture (n,7) and photo- 
dissociation (7,n) as the only reactions left in equilib- 
rium. 

This (n,"f) ^ equilibrium phase marks the on- 

set of the rprocess, as the balance between neutron cap- 
ture and photo-dissociation determines the equilibrium 
path of the rprocess. The rprocess may produce nuclei 
unstable towards fission at this time, whose fission prod- 
ucts may effect the r-process abundances. Once there are 
not enough free neutrons left to maintain (n, 7) (7, n) 
equilibrium, neutron capture and photo-dissociation re- 
actions freeze-out allowing the r-process nuclei to /3- 
decay back to stability, forming the progenitors to the 
observed r-process elements. 



problematic issues associated with the rprocess in the 
neutrino driven wind from the protoneutron star. Active- 
active neutrino transformation as well as active-sterile 
transformation through different channels has also been 
explored in Caldwell et al. 2000 42] and Nunokawa et 
al. 1997 [43. 

Here we consider matter enhanced transformation of 
neutrinos through the ^ Vs and Vf. ^ Vs channels 
as described by ^Jl- We define a sterile neutrino as 
one with interactions significantly weaker than the nor- 
mal weak interaction. These interactions must be weak 
enough as to not contribute significantly to the decay 
rate of the Z° boson Our results are independent 

of the model which produces a light sterile neutrino, as 
any SU(2) standard model singlet could be a light sterile 
neutrino. 

For the following equations we use for the electron 
neutrino wavefunctions and ^'s for the sterile neutrino 
wavefunctions. We consider evolution of flavor eigen- 
states in matter without the potential created by back- 
ground neutrinos in the form 45] 



III. MATTER-ENHANCED ACTIVE-STERILE 
NEUTRINO OSCILLATIONS 

McLaughlin et al. 1999 |^ demonstrated that active- 
sterile neutrino transformation could circumvent the 



7,^ 

or 



where 



0e(O VA 

\/A -0e(r) 



*e(r) 



(5) 



4 



) = ^{±2^GF[N-{r) - N+{r) - 



N„(r) 
2 



E 



-6m^ cos26'„) 



(6) 



and 



5m 



■ sin 29,, 



(7) 



In equation lO, the positive sign designates electron 
neutrinos and the negative sign electron anti-neutrinos. 
Here, Gf is the Fermi constant and N~ {r),N^ (r), and 
Nn{r) represent the number density for the electrons, 
positrons, and neutrons in the medium. The neutrino 
mixing angle is and the vacuum mass-squared split- 



ting term is defined as (5m 



where mi and 



7712 are the masses of the mass eigenstates. 

The potential is defined to be proportional to the net 
weak charge. 



V{r) = 2V2GF[N-{r) ~ N+{r) - 7V„(r)/2], 



meaning neutrinos of energy 



Eres{r) = ± 



^771^ cos 29y 



(8) 



(9) 



undergo a Mikheyev-Smirnov-Wolfenstein (MSW) reso- 
nance at a positive (i/g) or negative (pe) potential. Not- 
ing the electron fraction is 



N-{r)-N+{r) 
Np{r)^N4r) ' 



(10) 



we can immediately see how Yg influences the sign of the 
potential, 



(11) 



Electron neutrinos transform for a positive potential and 
electron anti-neutrinos for a negative potential. Above, 
p{r) is the matter density a distance r away from the 
protoneutron star, mN is nucleon mass, and Np(r) is the 
proton number density. Since the local material is elec- 



trically neutral we take Np 



N: 



N+. 



Using the above substitutions, we can rewrite equation 
© as 



Mr) = ± 



3GFp(r) 
2^/2mj\i 



Y,^-\-—cos29,. (12) 



We calculate equation ^ numerically for neutrino sur- 
vival probabilities passing through matter above the pro- 
toneutron star. We assume no flavor transformations oc- 
cur until the neutrinos leave the surface of the protoneu- 
tron star, hence all neutrinos reside in the active flavor 
eigenstates, v^. and V,,, with no sterile mixing, Us and 
Vg, at this time. We begin the neutrino transformation 
evolution a Fermi-Dirac neutrino and anti-neutrino spec- 
trum. As the neutrinos depart from the protoneutron 
star, the evolving survival probabilities from equation 
© cause the neutrinos to depart rapidly from this initial 
spectrum. We use electron neutrino and anti-neutrino 
luminosities of L^, = 10^^ ergs and Lp — 1.3 x 10^^ ergs 
respectively, electron neutrino and anti-neutrino temper- 
atures of = 3.5 MeV and Tp = 4.5 MeV respectively, 
and an effective chemical potential of zero. We chose 
these values to be representative of the general qualita- 
tive behavior of the expected neutrino flavor transforma- 
tion in this environment. For these luminosities, electron 
neutrinos and electron anti-neutrinos strongly determine 
the initial electron fraction, although additional influence 
on the electron fraction comes from electron and positron 
capture through the backward rates of the following re- 
actions: 



and 



(13) 



(14) 
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FIG. 2: 

sterile neutrinos above the surface of the neutron star by 
active-sterile MSW neutrino transformation. The active- 
sterile neutrino mixing parameters are 5m^ = 2eV'^ and 
s\v? 29v = 7 X 10~^, with the wind parameters of Fig. 
The neutrinos shown have energy iJi/ ~ 11 keV. 



Active-sterile MSW neutrino resonances occur when 
(j)e — 0, providing a direct relationship to Ye from the 
neutrino mixing parameters, 5w? and sin^ 29,^ and the 
electron neutrino and anti-neutrino energies. 



We now point out the main features of active-sterile 
neutrino mixing in this environment: the detailed behav- 
ior of Vf, ^ Vs mixing is discussed in [4lLl46|. As neutri- 
nos leave the surface of the protoneutron star, low energy 
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electron neutrinos are the first to undergo transformation 
to sterile neutrinos, Vg, due to the large potential, equa- 
tion 1)11(1 . Such a transformation is shown in Fig. [3 This 
reduces the efficacy on the forward rate of equation l|13|) 
lowering both the potential and the electron fraction, al- 
lowing progressively higher energy electron neutrinos to 
transform. Once the electron fraction drops below 1/3, 
the highest energy electron anti-neutrinos transform to 
sterile neutrinos due to the now negative potential. The 
potential continues to drop towards its minimum, after 
which it will ascend back towards its zero value. On its 
way back, upon the potential reaching this zero value, 
all of the electron anti-neutrinos above a threshold en- 
ergy convert from F^'s back to VeS. As the mass ele- 
ment reaches the alpha particle formation epoch, very 
few of the original electron neutrinos and most of the 
original electron anti- neutrinos are present. This severely 
impedes the forward rate of equation 1)12(1 due to the lack 
of z^e's, and also prevents the dramatic electron fraction 
rise associated with the a effect, allowing a neutron-rich 
neutron-to-proton ratio needed for r-process element for- 
mation. The suppression of the a effect is dependent on 
the neutrino mixing parameters, 5w? and sin^20^, and 
we describe later which parameters allow for r-process 
elements to be produced |41|. 



IV. NETWORK CALCULATIONS 

We calculate the nuclide abundances for a mass ele- 
ment leaving the surface of a protoneutron star by linking 
three network calculations: a nuclear statistical equilib- 
rium (NSE) code ^ full reaction rate network code 
|40| . and an r-process reaction network code 0, lislj . 
Throughout the calculation we self-consistently solve the 
neutrino evolution through the integration of equation 
(|5|). In all regimes we are guided by the structure of 
the protoneutron star from I49|, and from the wind as 
described in ^3 described in Sec. II. 

As the mass element first leaves the protoneutron star 
with TmcV ~ 3, the NSE code takes distance and density 
and calculates all relevant thermodynamic quantities in- 
cluding positron and electron number densities. We use 
NSE to calculate the neutron, proton, and heavy element 
[A > 40) abundances. We also compute the weak reac- 
tion rates, generating an updated electron fraction, Yg,. 
This allows us to self-consistently solve the matter en- 
hanced (MSW) equation and hence the neutrino survival 
probabilities. 

Once the mass element leaves the NSE regime, at 
TkeV ^ 900, we use the full reaction-rate network to 
calculate the evolution of the mass element. In this 
regime intermediate nuclide production is accounted for, 
A < 120, utilizing strong, electromagnetic, and weak in- 
teractions rates from |50j |. As before, we continue to use 
the updated electron fraction to self consistently calcu- 
late the MSW neutrino evolution. 

As the mass element reaches the neutron capture 



regime where charged particle reactions have frozen out, 
we use an r-process reaction rate code, for Tkev ^ 200. 
The relevant nuclear reactions included are neutron cap- 
ture, photo-disintegration, beta decay, charged-current 
neutrino interactions, and beta-delayed neutron emis- 
sion. We use the nuclear mass model from [sil. (3- 
decay rates from 52] and neutron capture rates from |50j |. 
While this mass model is known to produce a gap after 
the second peak, we choose it since it is the largest self- 
consistent data set available. We include the effects of 
nuclei undergoing fission as it is relevant in this regime. 



V. FISSION IN r-PROCESS 
NUCLEOSYNTHESIS 

It is important to understand which nuclides become 
unstable to fission and if and where fission terminates 
the r pro cess. Some previous works such as Cowan et al. 
1999 lUl and Meyer 2002 ^ have included termination 
of the rprocess through fission. Meyer 2002 '25| used a 
network termination boundary of Z = 91 from proton 
drip line to neutron drip line with only one isotope ^'^^Pa 
undergoing fission. Cowan et al. 1999 employed a 
probability of unity for spontaneous fission for nuclei with 
atomic weight greater than 256. Other calculations such 
as Goriely et al. 2004 js^l and Panov and Thielemann 
2004 [s^ , employed a set of fission rates over a range of 
heavy nuclei. 

Perhaps even more important than which nuclei un- 
dergo fission is the role the nuclide distribution of the 
daughter products play. There are two fission modes gen- 
erally considered, described by their daughter product 
distributions, symmetric and asymmetric. Symmetric fis- 
sion results in the progenitor nucleus being split nearly 
in half by fission, resulting in daughter products of the 
same nuclide for even numbers of protons and neutrons 
in the progenitor and adjacent nuclides for odd permu- 
tations of protons and neutrons. For asymmetric fission, 
the progenitor nucleus splits unevenly, resulting in one 
daughter product being proportionally heavier than the 
other. 

Cowan et al. 1999 [H and Rauscher et al. 1994 ^ 
employ symmetrical fission, to determine the daughter 
products, while Panov and Thielmann 2004 [s^l utilize 
a weighting function to determine the various daughter 
products. Ternary fission has been suggested as well by 
Ref. [sgI Isrf. At present, a consensus on the choice of 
fission daughter products has not been formed. 

In this work we go one step beyond a single nucleus ter- 
mination point. While we terminate the network path at 
A > 270 by taking the probability of spontaneous fission 
to be unity above this point, we also include beta-delayed 
fission as well for nuclei ranging from 81 < Z < 100 and 
140 < N < 164. Beta-delayed fission has the greatest 
importance during late stages of the rprocess [13|- Our 
approach is sufficient to illustrate the qualitative impact 
of ^ transform on the rprocess. For a detailed 
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comparison on a nucleus by nucleus basis with data, it 
would be preferable to include neutron-induced fission as 
well. In the future, such rates may become available [58l |. 

As stated above, the choice of fission daughter prod- 
ucts can have a significant impact on the distribution of 
final r-process element abundances. We examined two 
models, symmetric fission where the two daughter prod- 
ucts are taken to be as close to atomic weight and number 
as possible, and asymmetric fission using daughter prod- 
ucts suggested by Ref. |0, which takes the daughter 
products to be about 40% and 60% of the mass of the 
progenitor nucleus. 
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FIG. 3: Electron fraction, Ye is plotted vs. T^ev for the wind 
parameters in Fig. Q with active-sterile neutrino mixing pa- 
rameters of Sm^ — 10 eV^ and s'rn^ 20v = 0.001. The appro- 
priate neutrino interactions are turned off at Tkev = 850 for 
the neutrino-nucleon interactions only and no neutrino inter- 
actions cases. The continues to increase even after heavy 
element formation, due primarily to neutrino-nucleon inter- 
actions. 



VI. RESULTS AND DISCUSSION 

It is convenient to cast our results in the form of both 
the electron fraction, Yg = p/ {n+p), and abundance pat- 
terns to allow for direct comparison with previous theo- 
retical and experimental results. We begin with a calcu- 
lation of electron fraction, Y^, since this is a key indicator 
of r-process element production in the neutrino-driven 
wind environment and it facilitates comparison with pre- 
vious work over the neutrino mixing parameter region, 
Sm^ and sm^26^. We then discuss the implications of 
our results and compare our calculated abundance pat- 
terns with both solar data (52,i60l and halo star data |6ll| . 
We also compare neutrino mixing parameters favorable 
to the rprocess against the parameter region that will be 
probed by Mini-BooNE |63|. We pay particular atten- 
tion to the effect of fission cycling since this process is a 
determining factor in the behavior of the system. 
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FIG. 4: An r-process pattern (blue line) producing only 
the second and third of the three r-process peaks occurs in 
the neutrino-driven wind when there is active-sterile neutrino 
mixing. The dramatic impact of active-sterile neutrino oscil- 
lations on the rprocess is seen immediately when contrasted 
against the abundances produced without neutrino mixing 
(red line). For the blue line neutrino mixing parameters of 
Sm'^ = 2eV^ and sin^ 29^ = 7 x 10~^ are used together with 
same astrophysical conditions as in Fig. The general fea- 
tures of the r-process pattern above A > 130 found in solar 
system data [59] (black line) are reproduced when Ye < 0.3 
at the onset of neutron capture element formation. The so- 
lar data is scaled to the simulation such that the sum of the 
A = 195 nuchdes is 2.3 x 10"^. 



A. Variation of Neutrino Mixing Parameters 

A key indicator of whether the rprocess will occur is 
the electron fraction, Yf,, at the onset of neutron cap- 
ture element formation. As described earlier, the a ef- 
fect together with low entropy is a prominent cause of the 
reduction of free neutrons, stifling the rprocess. Active- 
sterile neutrino oscillations allow the a effect to be cir- 
cumvented during the epoch of a particle formation as 
electron neutrinos are converted to their sterile counter- 
parts, preventing electron neutrino capture on neutrons, 
and producing a neutron-rich environment. 

Our reaction rate network self-consistently accounts for 
the evolution of Yg throughout the duration of nucleosyn- 
thesis of the mass element. We begin by comparing Y^. 
with a previous calculation. Fig. 7 of Ref. [41|, at the 
onset of heavy element formation (Fig. QJd) and find that 
good agreement is reached. Then we further track the 
Ye up to the point of neutron capture element formation, 
TkeV ~ 200, and we find that for Y^ < 0.35 at this time, 
the environment is sufficiently neutron-rich for a rprocess 
to occur (Fig. ^). 

We note there is a slight increase of Yg between the two 
times, as shown between the two panels of Fig. ^ Since 
charged-current reactions affect the Yg as the mass el- 
ement moves farther away from the protoneutron star. 
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FIG. 5: The general features of the r-process pattern found in 
halo stars are reproduced in the neutrino-driven wind when 
active-sterile neutrino oscillations occur. Only the second, 
A 130 {Z « 55), third peaks, A ki 195 {Z ^ 80), as well 
as the rare earth bump inbetween, are reproduced. The halo 
star CS 22892-052 |6(l | abundances are plotted for compar- 
ative purposes (dots) along with the solar system r-process 
abundances from |5fll ] (black line). Neutrino mixing param- 
eters of 6m^ = 2eV^ and sin^ 29v — 7 x 10~^ were used for 
the neutrino driven wind abundances (blue line) , although the 
pattern shown is fairly insensitive to the exact choice of mix- 
ing parameter (see Fig. [SJ. The astrophysical conditions are 
the same as in FigQ The measured abundances are scaled to 
the simulation such that the sum of the Z = 70 isotopes in 
each abundance curve have an abundance of 10^'^. 

the increase in Yg must be the result of one or more 
charged-current reactions, electron neutrino capture on 
neutrons (the forward rate of equation l|13() ) increases 
Ye, making the environment less neutron-rich. Although 
electron neutrinos that undergo transformation to sterile 
neutrinos do not affect the Y^, not all electron neutri- 
nos undergo transformation. The effect of these residual 
electron neutrinos is shown in Fig. 13 where we plot 
the electron fraction between heavy nucleus production 
and neutron capture element formation, for the neutrino 
mixing parameters of ^m^ = 10 eV^ and sin^ 29^ = 0.001 
Electron neutrino capture on neutrons is therefore re- 
sponsible for the increase in seen in Fig. ^ 



B. i?-process Abundance Patterns 

In Fig. ^ we examine the consequences active-sterile 
neutrino oscillations have on abundance patterns in the 
neutrino-driven wind environment and find conditions 
suitable for the rprocess in this environment for Ye ^ 
0.35. We note that none of our calculations reproduce the 
first r-process peak at A ~ 80. The resulting abundance 
pattern features the second and third r-process peaks and 
we note that the deficit after the second r-process peak, 
A ~ 130 is a consequence of the mass model we are using. 



Fig. 0] (blue line). The solar abundance of A = 195 is 
scaled to 2.3 x 10""^ for comparison. 

We classify our abundance patterns by Y, at the onset 
of neutron capture element formation '63'|. For 0.35 > 
Ye > 0.30, there are too few free neutrons to produce the 
third r-process peak, A ^ 195, with the same proportions 
as the second r-process peak, A ~ 130, as is indicated by 
observational data. Electron fractions of Ye < 0.3 lead 
to conditions sufficiently neutron-rich to have the second 
and third r-process peaks in similar proportions to the 
observational data. Ye lower than Ye < 0.25 produces 
r-process elements around A > 250 that are strongly 
susceptible to fission, which can in turn repopulate re- 
gions at and around the ^ ~ 130 peak. We will return 
to this fission cycling in the following section. Fig. 
shows our abundance pattern plotted as a function of 
proton number, further demonstrating the second and 
third r-process peak production. The three abundances 
are scaled such that Z = 70 has an abundance of 10^"^. 

Early calculations reproduced the general features of 
the solar system r-process pattern by constructing a final 
r-process abundance pattern through a mass weighted 
average of individual abundances patterns resulting mul- 
tiple mass trajectories d, Q, I3- have investigated 
different wind conditions, including changes in neutrino 
luminosity and outflow timescale. The effect of varying 
the outflow timescale is shifting the range of neutrino 
mixing parameters, Sm^ and sin^ 20, for which a success- 
ful rprocess occurs. The effect of reducing the neutrino 
luminosity is contracting the range of the neutrino mix- 
ing parameters yielding a successful rprocess. While we 
find the absence of an ^ ~ 80 peak for a wide range 
of conditions, it would be necessary to calculate a time 
integrated history over a particular series of wind condi- 
tions to determine in a given neutrino driven wind model 
produces a significant first r-process peak. 



C. Effects of fission cycling 

The choice of neutrino mixing parameters affects the 
free neutron abundance and hence the electron fraction 
at the time of neutron capture element formation. We 
find that fission cycling is a mechanism for generating a 
stable r-process pattern once a threshold Ye is reached. 
Ye ^ 0.25 in our environment. All neutrino mixing pa- 
rameters which are below this Ye produce the same gen- 
eral r-process pattern, shown in Fig. 

Fission cycling behaves as a mechanism linking the sec- 
ond and third peaks of the rprocess. As very heavy nu- 
clides are produced, A > 250, they may become unstable 
towards fission and provide an effective termination point 
to the r-process path. The resulting daughter nuclides 
from this fission repopulate abundances of the second r- 
process peak. 

The consistency of the r-process pattern is dependent 
on the equilibrium path of the rprocess. Early studies 
analyzing the mass numbers of both the peaks and the 
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FIG. 6: Fission cycling produces a consistent overall r-process pattern in the neutrino-driven wind, as long as Ye < 0.3 at 
the onset of neutron capture element formation (right panel). In the left panel we show neutron abundance vs. TkeV for three 
sets of neutrino mixing parameters, 5m^ = 10 eV^ and sin^ 26v = 10~^, Sm^ — 5eV^ and sin^ 26^1, — 10~^, and 5m^ = 2eV^ 
and sin^ 26'„ = 7 x 10~^ corresponding to electron fractions of = 0.1174, Ye = 0.1752, and = 0.2498 respectively. The 
Ve in the legend occurs at the start of heavy element formation, Tkev ~ 200. Lower Ve's lead to more neutrons but a similar 
abundance pattern. In the right panel we show the abundance patterns from the same mixing parameters as used in the left 
panel. The abundances curves and astrophysical conditions are the same as figure 



closed neutron shells, yielded a r-process path well away 
from the neutron-drip line. In general the equilibrium 
path of the rprocess is determined by the balance be- 
tween neutron capture and photo-dissociation for partic- 
ular temperature and density conditions. In our calcula- 
tions the environment becomes sufficiently neutron-rich 
to generate a r-process path along the neutron-drip line, 
resulting in a consistent r-process abundance pattern of 
the second and third peaks. 

In the remainder of this section we discuss some details 
about the effects of fission. It is difficult to make any 
definitive comments on the specific effects that sponta- 
neous fission and beta-delayed fission have on particular 
final r-process nuclide abundances due to uncertainties in 
fission models; however, some general comments are war- 
ranted. The impact of fission is greatest in the regimes 
around the A ~ 130 peak and past the A ^ 195 peak, the 
regimes where fission deposition and depletion occurs. 

The details of the pattern around the A ~ 130 peak 
is caused by the deposition of fission daughter products, 
which in turn are strongly affected by the choice of fission 
daughter product distribution, in the region of 120 < A < 
140. Two modes of daughter product distributions that 
we examined are symmetric and asymmetric as described 
in Ref. and ^^l- Symmetric fission modes generally 
produce an abundance pattern similar to that produced 
by an r-process calculation without fission, since material 
is mostly deposited near the closed neutron shell region of 
the A ^ 130 peak; however, asymmetric modes can make 
significant material deposits into the A 100 region and 



erase the deficit past the A ~ 130 peak. It has been 
suggested that fission could form elements yl < 120 in 
metal poor stars previously, e.g. Ref. [63 |. 

Abundance structure of the depletion region past the 
A ~ 195 peak is strongly affected by several factors. 
These factors are the amount of material that reaches 
this region, how strongly fission depletion occurs, and for 
what elements it occurs. Most previous work on fission in 
the rprocess involves calculations which effectively termi- 
nate the rprocess using a 100% fission probablity in the 
regime of 260 S ^ ^ 280. Given that the nuclear physics 
of such neutron-rich nuclei is not well understood, it is 
difficult to draw conclusions about any detailed features 
in this region. Nevertheless, fission cycling produces a 
distinct and stable abundance pattern with overall fea- 
tures that are in agreement with data. 

D. Comparison with Halo Star Abundances 

Halo stars, by virtue of their ancient nature, give us 
a glimpse of the abundances from the earliest r-process 
progenitors. The similarity between the r-process pat- 
tern in halo stars and that seen in the solar rprocess sug- 
gests a single extremely consistent rprocess mechanism is 
operating for the heavier r-process elements. We present 
our abundance pattern for comparison with the halo star 
CS 22892-052 59, 61] and the solar system rprocess H^l 
in Fig. Inland Fig. |^ There is rough agreement between 
our data. Fig. (black line), and CS 22982-052, Fig. 
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0.001 , 0.1 

sin'(2e^) 

FIG. 7: Contours indicating the parameter space for which 
we find an rprocess (defined as Ye < 0.35) in the prospective 
neutrino mixing parameter range of Mini-BooNE. The light 
gray region signifies an rprocess for an outflow timescale of 
r = 0.1s, and similarly the dark gray region for r — 0.3 s, 
with a reminder of the astrophysical conditions listed in Fig. 
El The blue region is the combined NSBL and LSND 90% CL 
allowed region [62l . The parameter space which produces a 
successful rprocess is much larger than the space probed by 
experiment. 



(dots), between 50 < Z < 80 the region of the second 
and third r-process peaks. The gap after the Z ~ 50 
peak is the same as the gap after A ~ 135 in the pre- 
vious plot, shown in Fig. |S1 and is feature of the mass 
model we are using. Our final abundances are allowed to 
decay over a timescale of 14 Gyr for consistency with the 
accepted age of CS 22982-052, using decay modes and 
lifetimes from Ref. [s^. Adjusting the decay timescale 
primarily affects the abundances of Th^^^ and U^'^* since 
their half- lives are the the shortest of the long lived quasi- 
stable elements, nearly the same order as the halo star. 
The remaining quasi-stable elements generally have de- 
cay paths to nearby stable nuclides so their decay does 
not produce a noticeable effect to the abundance pattern. 



E. Comparison with Combined NSBL and LSND 
Data 

The Mini-BooNE experiment is designed to study neu- 
trino mixing in the range of the positive LSND result. 
One possible implication of this result is the existence of 
a sterile neutrino. In such a scenario the range of mixing 
parameters yielding a successful rprocess in our calcula- 
tions overlaps with the parameter range of Mini-BooNE. 
In Fig. [71 we compare our r-process allowed region for 
two outflow timescales, r = 0.1s and r — 0.3 s, with 
the combined null-short baseline (NSBL) and liquid scin- 
tillator neutron detector (LSND) (3-1-1) 99% confidence 
level (CL) data, with lower CL regions naturally resid- 



ing inside the 99% region. A compilation of the NSBL 
and LSND data is described in detail in Ref. ^2]. It 
should be noted that several factors impact the location 
of the r-process region and hence any compatibility with 
the NSBL and LSND data. The choice of expansion 
timescale, t, alters the electron fraction and shifts the 
points where active-sterile MSW flavor transformation 
occurs. 

Finally, we caution that the inclusion of the neutrino 
background terms in the weak potential can significantly 
change the r-process abundances ^^|6^|63. However, 
at late times the luminosity decreases, decreasing the ef- 
fect of background terms [46L l6^ . In our calculations 
which do not include the background, we find that a re- 
duced electron neutrino luminosity oi ~ 5 x 10^° ergs 
and a reduced electron anti-neutrino luminosity of Lp — 
6.5 X 10^° ergs can still yield a robust r-process solution. 
Future studies of neutrino transformation fully inclusive 
of background terms done together a hydrodynamical 
wind calculation will determine the extent to which r- 
process nucleosynthesis products will be produced and 
also escape the star. 



VII. CONCLUSIONS 

We have examined nucleosynthesis in the neutrino- 
driven wind in the post-core bounce supernovae envi- 
ronment. We find that inclusion of active-sterile neu- 
trino transformation yields an r-process pattern where 
only the second and third of the r-process peaks are pro- 
duced for a broad range of neutrino mixing parameters. 
Analysis of data from halo stars has suggested that there 
are at least two separate types of events which make r- 
process nuclei 01. Further- more the analysis suggests 
one of these mechanisms operates early in the evolution 
of the universe and makes primarily the second and third 
r-process peaks. 

Fission cycling is responsible for the establishment of a 
stable r-process pattern. As long as many neutrons exist 
during the phase of neutron capture, nuclei are created 
with sufficiently high mass number that they are unstable 
to fission. Below a threshold (Ye ^ 0.3) in our model 
the resulting abundance pattern is not sensitive to the 
magnitude of the free neutron abundance and therefore 
the number of times the fission cycle occurs. 

We find our allowed r-process region has overlap with 
the parameter space to be probed by the Mini-BooNE 
experiment. The extent of this overlap is dependent on 
uncertainties in the astrophysical environment and in the 
neutrino transformation calculations; however, the exis- 
tence of a sterile neutrino in the parameters space that 
will be probed by Mini-BooNE will directly influence el- 
ement production in core-collapse supernova. 
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